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Overview 


□ Supersonic Commercial Transport Aircraft Design 

□ Game Change Approach in Aircraft Design 

□ Multidisciplinary Design Optimization tool 

□ Multidisciplinary Analysis of the Baseline Configuration 
♦> Structural and Aerodynamic Models 

♦> Modal Analyses 

❖ Flutter Analyses 
♦> Trim Analyses 

♦> Landing and Ground Control Loads 

❖ Buckling and Strength Analyses 

□ First Optimization Run 

□ Second Optimization Run 

□ Third Optimization Run 

□ Conclusions 

□ Future Studies 
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Supersonic Commercial Transport Aircraft Design 
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Gulf Stream 


□ 


❖ 


❖ 


Major Issues 

Safety 

> Light weight airframe can cause strength, buckling, aeroelastic, and aeroservoelastic problems. 
Sonic boom 

> Supersonic flight of "commercial transport" aircraft allowed only over the ocean. 

> Perceived Loudness in decibels 
S NASA's N+2 goal: 85 PLdB 
S Concorde: 104 PLdB 

•S High Speed Civil Transport (HSCT): 99 PLdB 
Fuel efficiency 

> Light weight airframe 

> Reduced drag 

Developing N+2 Low-boom Supersonic Commercial Transport (LSCT) aircraft 

^ oe * n S Lockheed Martin 

Lockheed Martin: 79 PLdB 


❖ 


HSCT 


❖ 


❖ 


Gulf Stream 


❖ Aerion with "Airbus" 








Concorde 


Game Change Approach in Aircraft Design 


Problem Statement 


□ Design innovations are needed to further down the weight of an aircraft which current design technologies can take care of. 

Long Term Objective 

□ Use aeroelastic tailoring theory and active flexible motion control technique to satisfy the overall strain, aeroelastic, and 
aeroservoelastic instability requirements within given flight envelopes 

□ Use curvilinear sparib concept as well as composite ply angles for aeroelastic tailoring 


A pproach 


□ Simultaneously update structural as well as control design variables during early design phase 



❖ Perform topology optimization with curvilinear sparibs 

❖ Use aeroelastic tailoring up to V L 

❖ Use aeroservoelastic tailoring between V L and 1.15 V L 

Current Study 

□ Optimize baseline aircraft model 

❖ Use Lockheed Martin's configuration 

❖ Use aeroelastic tailoring up to 1.15 V L 

Curvilinear sparibs 



Mach Number 
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Multidisciplinary Design Optimization tool 


□ Based on Object-Oriented Optimization tool 

Open MDAO, Model Center, Visual Doc, etc. 


❖ 


Update 


O : Incorporated 
O : Modules are being developed 
# : Not included yet 


MSC/NASTRAN sol 103 
Total weight, CG location, 
mass moment of inertia, 
Frequencies, & mode shapes 


MSC/NASTRAN for small weight 
In-house code for large weight 

PI W = W T 


MSC/NASTRAN sol 105 
In-house code for computing MS 
Use safety factor of 1.5 


MS = 


Failure Load 


Design Load x Safety Factor 
PI S = -min(MS) 


In-house code for computing BLF 
Use safety factor of 1.5 


PI B = (1/2) 2 — (positive min(BLF) — 1/2} 



ZAERO code for flutter analyses 
In-house code for flutter speed 
tracking Vp 

PI F = 1 - 


Speed V 
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In-house code for 
massaging splined loads 
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Multidisciplinary Analysis of the 
Baseline Configuration 


Structural and Aerodynamic Models 



□ MSC/NASTRAN structural model 
❖ Total number of grids: 55,635 



□ ZAERO unsteady aerodynamic model 

❖ 5,060 surface elements 

❖ Six Mach numbers: 0.66, 0.89, 1.41, 1.80, 2.00, and 2.30 

❖ Sixteen reduced frequencies: 0., 0.005, 0.01, 0.05, 0.10, 0.15, 0.20, 
0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.65, 0.80, 1.0 
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Modal Analyses of the Baseline Configuration 


□ Based on six configurations 
♦> Gear up 

> DTOW(Design Take Off Weight), FFEP(Full Fuel Empty Payload), M2W(Mach 2 Weight), & ZFW(Zero Fuel Weight) 

S DTOW=FFFP (Full Fuel Full Payload) 

•S ZFW=EFFP (Empty Fuel Full Payload) 

♦♦♦ Gear down 

> DTOW(or FFFP) and DLW(Design Landing Weight) 


Mode 

Number 

Natural Frequency (Hz) 

Notes 

Gear-up 

Gear-down 

DTOW 

FFEP 

M2W 

ZFW 

DTOW 

DLW 

7 

2.049 

2.055 

2.071 

2.266 

2.048 

2.158 

Aft fuselage torsion 

8 

2.235 

2.262 

2.277 

2.554 

2.238 

2.424 

First symmetric fuselage bending 

9 

2.498 

2.509 

2.539 

2.993 

2.503 

2.714 

First symmetric wing bending 

10 

2.754 

2.769 

2.935 

3.415 

2.752 

3.265 

First anti- symmetric wing bending 

11 

3.060 

3.069 

3.115 

3.731 

3.057 

3.403 

Symmetric tail bending 

12 

3.562 

3.608 

3.689 

4.044 

3.574 

3.945 

Forward fuselage lateral bending 

13 

4.440 

4.449 

4.511 

4.790 

4.429 

4.602 

First anti- symmetric tail bending 

14 

4.456 

4.537 

4.555 

5.532 

4.437 

5.142 

Second symmetric wing bending 

15 

4.818 

4.842 

5.146 

5.832 

4.809 

5.542 

Second anti-symmetric wing bending 

16 

5.449 

5.465 

5.550 

6.158 

5.444 

5.994 

Symmetric aft inner wing bending 
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Flexible Mode Shapes (gear up: DTOW) 



Aft fuselage torsion 
Z nodal lin 



First symmetric fuselage bending 



First symmetric wing bending 
Z nodal line 


1 st Mode: 2.049 Hz 


2 nd Mode: 2.235 Hz 


3 rd Mode: 2.498 Hz 


First anti-symmetric wing bending Symmetric tail bending Forward fuselage lateral bending 



, n 4 th Mode: 2.754 Hz 

Structural Dynamics Group 


5 th Mode: 3.060 Hz 


6 th Mode: 3.562 Hz 
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Flutter Analyses of the Baseline Configuration 


❖ 


□ Based on four structural configurations at six Mach numbers 
Gear up 

> DTOW, FFEP, EFEP, & ZFW 

> Mach numbers of 0.66, 0.89, 1.41, 1.80, 2.00, & 2.30- 



O : 

O: 

1100 

1000 

900 

800 

? 700 

<D 

* 600 
■§ 500 
& 400 
300 
200 
100 
0 


Empty Fuel Empty Payload (EFEP) A : Full Fuel Empty Payload (FFEP) 
Zero Fuel Weight (ZLW) □: Design Take-Off Weight (DTOW) 

-10k S.L. 10k 20k 


■ 

1 1 1 1 1 

Based on 1% 
damping 
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Trim Analyses of the Baseline Configuration 

□ Trim flight conditions 


Load 

Case 

ID 


100 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


1100 


1200 


1300 


1400 


1500 


1600 


1700 


1800 


2100 


2200 

2300 

2400 


Maneuver 


Pull up 


Push over 


Pull up 


Pull up 


Push over 


Pull up 


Pull up 


Push over 


Pull up 


Push over 


Steady roll 


Abrupt roll 


Steady roll 


Abrupt roll 


Landing 


Cruise 


Gust Loads 


Landing 


Steady roll 


Abrupt roll 
Pull up 
Pull up 


Load 

Factor 


2.5g 


lg 


2.5g 


2.5g 


■lg 


2-5g 


2-5g 


■lg 


2.5g 


■lg 


0g 


0g 


L67g 


L67g 


lg 


lg 


2-7g 


lg 


0g 


0g 

0g 

1.67g 


Mach 

Number 


0.66 


0.66 


0.48 


2.00 


2.00 


1.41 


0.66 


0.66 


2.00 


2.00 


0.48 


0.48 


0.48 


0.48 


0.3092 


1.80 


0.89 


0.3092 


0.48 


0.48 

0.48 

0.48 


Weight 


DTOW 


DTOW 


DTOW 


M2W 


M2W 


DTOW 


ZFW 


ZFW 


ZFW 


ZFW 


DTOW 


DTOW 


DTOW 


DTOW 


DTOW 


DTOW 


ZFW 


DLW 


DTOW 


DTOW 

DTOW 

DTOW 


Landing 

Gear 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Up 


Down 


Up 


Up 


Down 


Up 


Up. 

Up. 

Up 


Altitude 


SL 


SL 


SL 


49,770ft 


49,770ft 


49,770ft 


SL 


SL 


49,770ft 


49,770ft 


SL 


SL 


SL 


SL 


SL 


55,000ft 


20,000ft 


SL 


SL 


SL 

SL 

SL 


Trim Variables 


BF(R=L) 


BF(R=L) 


BF(R=L) 


BF=TEF(R=L) 


BF(R=L) 


BF=TEF= AIL 1 = AIL2(R=L) 


BF(R=L) 


BF(R=L) 


BF=TEF(R=L) 


BF(R=L) 


Load Case 2100+2300 


Load Case 2200+2300 


Load Case 2100+2400 


Load Case 2200+240 1 


BF(R=L) 


BF=TEF(R=L) 


BF(R=L) 


BF(R=L) 


AIL1— AIL2(R— L) Aileron#2(L) 


AIL1=AIL2(R=-L) 

BF(R=L) 

BF(R=L) 
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Trim Analyses of the Baseline Configuration (continue) 


□ Trim results 

❖ In general trim angles are larger than NASTRAN 
results. 




Load Case 

100 

200 

300 

400 

500 

600 

700 

800 

900 

Trim Analysis 

Symmetric 

Nx (G) 

-0.007 

-0.003 

-0.005 

0.001 

0.002 

-0.004 

-0.016 

-0.006 

0.00; 

Nz (G) 

2.5 

-1.0 

2.5 

2.5 

-1.0 

2.5 

2.5 

-1.0 

2.5 

Pdot (rad/s 2 /g) 

None 

None 

None 

None 

None 

None 

None 

None 

None 

Qdot (rad/s 2 /g) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Pb/2V (rad) 

None 

None 

None 

None 

None 

None 

None 

None 

None 

Qc/2V (rad) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

a(°) 

7.75 

-2.50 

14.37 

8.32 

-2.81 

16.90 

4.10 

-1.04 

5.07 

Body Flap (°) 

2.01 

-6.07 

6.12 

-5.25 

5.42 

-25.68 

-8.01 

-1.59 

-12.9: 

Trailing-Edge Flap (°) 




-5.25 


-25.68 



-12.9: 

Aileron #1 (°) 






-25.68 




Aileron #2 (°) 






-25.68 




Mach Number 

0.66 

0.66 

0.48 

2.00 

2.00 

1.41 

0.66 

0.66 

2.00 

Altitude (ft) 

SL 

SL 

SL 

49770 

49770 

49770 

SL 

SL 

49771 

Weight Configuration 

DTOW 

DTOW 

DTOW 

M2W 

M2W 

DTOW 

ZFW 

ZFW 

ZFW 

Gear Configuration 

Up 

Up 

Up 

Up 

Up 

Up 

Up 

Up 

Up 

Load Case 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

Trim Analysis 

Sym. 

Asymmetric (sym. + Anti-sym.) 

Symmetric 

Nx (G) 

0.004 

-0.002 

-0.002 

-0.004 

-0.004 

-0.002 

0.002 

-0.022 

-0.00 

Nz (G) 

-1.0 

0.0 

0.0 

1.67 

1.67 

1.0 

1.0 

2.7 

1.0 

Pdot (rad/s 2 /g) 

None 

0.0 

0.0014 

0.0 

0.0014 

None 

None 

None 

None 

Qdot (rad/s 2 /g) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Pb/2V (rad) 

None 

0.0410 

0.0 

0.0410 

0.0 

None 

None 

None 

None 

Qc/2V (rad) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

a(°) 

-1.44 

0.40 

0.40 

9.74 

9.74 

13.91 

6.02 

4.95 

9.07 

Body Flap (°) 

11.89 

-2.86 

-2.86 

3.14 

3.14 

8.00 

-9.87 

-13.20 

22.5^ 

Trailing-Edge Flap (°) 







-9.87 



Aileron #1 (°) 


19.07 

48.63 

19.07 

48.63 





Aileron #2 (°) 


19.07 

48.63 

19.07 

48.63 





Mach Number 

2.00 

0.48 

0.48 

0.48 

0.48 

0.3092 

1.80 

0.89 

0.309 

Altitude (ft) 

49770 

SL 

SL 

SL 

SL 

SL 

55000 

20000 

SL 

Weight Configuration 

ZFW 

DTOW 

DTOW 

DTOW 

DTOW 

DTOW 

DTOW 

ZFW 

DLW 

Gear Configuration 

Up 

Up 

Up 

Up 

Up 

Down 

Up 

Up 

Dowi 
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Landing and Ground Control Loads of the 
Baseline Configuration 


□ 

□ 


Landing loads 
Ground control loads 


Ground Control Loads 


□ Emergency landing loads (applied to three engine structures] 

9G forward loading; 1.5G rearward loading; 3G sideway loading; & 
6G downward loading Landing Loads 


❖ 


Load Cases 


Case Number 


Load 


Right-MLG 


Left-MLG 


NLG 


Load 

Cases 


Static 

Condition 


3-point 
braked roll 


2-point 
braked roll 


Case Number 


3009 (DTOW) 
& 4009 (DLW) 


3010 (DTOW) 
& 4010 (DLW) 


Load 


FX 


FY 


FZ 


FX 


FY 


FZ 


FX 


FY 


FZ 


Right-MLG 


0.0 


0.0 


FM 5 , 

^•^^CG2NG^T 


X NG2MG 


0.8FM 


st 


0.0 


FM 


st 


0.8FM 


st 


0.0 


FM 


£L_ 


Left-MLG 


0.0 


0.0 


FM 


st 


0.8FM 


st 


0.0 


FM 


st 


0.8FM 


st 


0.0 


FM 


£L_ 


NLG 


0.0 


0.0 


FN 5t = 


X CG2MG 


W T 


X NG2MG 


0.0 


0.0 


2d 


CG2MG 


FM 5t + 2 (Z CG — Z MGCP ) 0.8FN 


X GG2NG 


0.0 


0.0 


0.0 


Level + 
Trim load 


3001 (DTOW) 
& 4001 (DLW) 


FX 


0.25FM 


Lv 


0.25FM 


Lv 


0.25FN 


Lv 


FX 


0.8FM 


st 


0.8FM 


st 


0.0 


FY 


0.0 


0.0 


0.0 


FY 


0.0 


0.0 


0.0 


FZ 


FM^ = f TMn W, 


LMG ' 


FM 


Lv 


FN i£ = UW, 


LNG \ 


Spin up + 
Trim load 


3002 (DTOW) 
& 4002 (DLW) 


FX 


(0.8 x 0.8)FM 


Lv 


(0.8 x 0.8)FM 


Lv 


(0.8 x 0.8)FN 


Lv 


FY 


0.0 


0.0 


0.0 


Dynamic 

roll 

braking 


W r 


3011 (DTOW) 
&4011 (DLW) 


, J CG2MG ■ . _s_ ..u 

a NG2MG L a NG2MG + 

where, E = {. Z CG - Z NGCP - (X rG - 

Xngcp)$} 

k Z MGCP ~ Z NGCP 


+ ' 


fd 


CG2NG 


liE 


FZ 


FM 


st 


FM 


st 


FZ 


0.8FM 


Lv 


0.8FM 


Lv 


0.8FN 


Lv 


Spring back 
+ Trim load 


3003 (DTOW) 
& 4003 (DLW) 


FX 


-(0.8 X 0.8)FM 


Lv 


-(0.8 X 0.8)FM 


Lv 


-(0.8 X 0.8)FN 


Lv 


and 5 — ( X CG — X NGCP ) 


Kmgcp -Xngcp 


FY 


0.0 


0.0 


0.0 


FZ 


0.8FM 


Lv 


0.8FM 


Lv 


0.8FN 


Lv 


Lateral 
drift + 
Trim load 


3004 (DTOW) 
& 4004 (DLW) 


FX 


(0.4 x 0.75)FM,,„ 


(0.4 x 0.75)FM 


Lv 


0.4FN 


Turning 

Condition 


3012 (DTOW) 
& 4012 (DLW) 


FX 


0.0 


0.0 


0.25FN 


st 


FY 


0.5FM 


SL. 


0.5FM 


SL. 


0.5FN 


SL. 


Lv 


FZ 


FM 


st_ 


FM 


st_ 


FN 


5L_ 


FY 


(0.25 x 0.75)FM 


Lv 


(0.25 x 0.75)FM 


Lv 


0.25FN 


Lv 


FZ 


0.75FM 


Lv 


0.75FM 


Lv 


FN 


Lv 


Right one 
gear + Trim 
load 


3005 (DTOW) 
& 4005 (DLW) 


FX 


0.25FM 


Lv 


0.0 


0.0 


Nose wheel 
yaw & 
steering (1) 


3013 (DTOW) 
& 4013 (DLW) 


FX 


0.0 


0.0 


0.0 


FY 


0.0 


0.0 


0.8FN 


£L= 


FZ 


FM 


£L_ 


FM 


£L_ 


FN 


5L_ 


FY 


0.0 


0.0 


0.0 


FZ 


FM 


Lv 


0.0 


0.0 


Left one 
gear + Trim 
load 


3006 (DTOW) 
& 4006 (DLW) 


FX 


0.0 


0.25FM 


Lv 


0.0 


Nose wheel 
yaw & 
steering (2) 


FX 


0.8FM 


£L_ 


0.0 


0.0 


3014 (DTOW) 
& 4014 (DLW) 


FY 


0.0 


0.0 


0.0 


FZ 


FM 


st 


FM 


2dcG2M£FM 5t + ( Z CG — Z MGCP ) 0.8FM 


st 


FY 


0.0 


0.0 


0.0 


X CG2NG 


FZ 


0.0 


FM 


Lv 


0.0 


Side load 
RtoL + 
Trim load 


3007 (DTOW) 
& 4007 (DLW) 


FX 


0.0 


0.0 


0.0 


FY 


(0.8 x 0.5)FM,,„ 


(0.6 x 0.5)FM,,„ 


0.0 


Nose wheel 
yaw & 
steering (3) 


FX 


0.0 


0.8FM 


SL. 


0.0 


3015 (DTOW) 
& 4015 (DLW) 


FY 


0.0 


0.0 


0.0 


FZ 


FM 


st 


FM 


2dcc2MC FMst + (Zee — Z^ccp)0-8FM 


st 


FZ 


0.5FM 


Lv 


0.5FM 


Lv 


0.0 


±CG2NG 


Side load 
LtoR + 
Trim load 


3008 (DTOW) 
& 4008 (DLW) 


FX 


0.0 


0.0 


0.0 


FY 


-(0.6 x 0.5)FM ; , 17 


-(0.8 x 0.5)FM ; , 1; 


0.0 


Reversed 

braking 


3016 (DTOW) 
& 4016 (DLW) 


FX 


-0.55FM 


st 


-0.55FM 


st 


0.0 


FY 


0.0 


0.0 


0.0 


FZ 


FMc 


FMc 


0.0 


FZ 


0.5FM 


Lv 


0.5FM 


Lv 


0.0 


DTOW 


fLMG-0-36; ^lng —0*0639; Trim load case ID = 1500 


2G Taxi 


DLW 


fLMG-1-20; fuvG=0.1477; Trim load case ID = 1800 


3017 (DTOW) 
& 4017 (DLW) 


FX 


0.0 


0.0 


0.0 


FY 


0.0 


0.0 


0.0 


FZ 


2FM 


st 


2FM 


st 


2FN 


st 
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Buckling and Strength Analyses 


□ Based on five analysis sets 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Load Cases 

1 

Up 

DTOW 

100, 200, 300, 600, 1100, 1200, 1300, 1400, & 1600 

2 

Up 

ZFW 

700, 800, 900, 1000, & 1700 Center Line 

3 

Up 

M2W 

400 & 500 ^ 

4 

Down 

DTOW 

3001 ~ 3017 + 3018 ~ 3021 (emergency) + 1500 (for 
landing) 

5 

Down 

DLW 

4001 -4017 + 4018 ~ 4021 (emergency) + 1800 (for 
landing) 


□ Minimum buckling load factors from each analysis set 



Wing 



Case 300; DTOW 


El 


Passenger Floor 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Case 

Number 

Load Case 

Minimum 
Buckling Load 
Factor 

Buckling 

1 

Up 

DTOW 

300 

2.5G pull up; M=0.48 

0.152 

I ye&Bi 

2 

Up 

ZFW 

1700 

2.7G gust loads; M=0.89 

0.195 

1 y es l 

3 

Up 

M2W 

400 

2.5G pull up; M=2.00 

0.151 

\vey 

4 

Down 

DTOW 

3006 

Left one gear landing 

1.71 

no 

5 

Down 

DLW 

4006 

Left one gear landing 

1.52 

no 



Main Landing Gear Bay 


Center Engine 


Buckling Load Factor > 1 or 

Buckling Load Factor < 0 : requirement 


□ Minimum margins of safety from each analysis set 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Case 

Number 

Load Case 

Minimum 
Margin of 
Safety 

Failur 

1 

Up 

DTOW 

1400 

1.67G abrupt roll; M=0.48 

-0.999 

/ ye^V 

2 

Up 

ZFW 

1700 

2.7G gust loads; M=0.89 

-0.998 

/ Y es \ 

3 

Up 

M2W 

400 

2.5G pull up; M=2.00 

-0.997 

I yes 1 

4 

Down 

DTOW 

3013 

Nose wheel yaw & steering (1) 

-0.781 

\ y es / 

5 

Down 

DLW 

4003 

Spring back landing 

-0.657 

\ye/ 


Margin of safety > 0 : requirement 
Safety factor =1.5 

Failure Load 

MS = 

Design Load x Safety Factor 
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First Optimization Run 




Design Variables for the First Optimization Run 


□ 


Baseline configuration is in infeasible domain 
Manually increase 111 ply thicknesses 

> Wing, tail, inner-wing, and fuselage skins 
S Three variables 

• 1 st , 3 rd ,& 2 nd = 4 th 

> Wing, inner-wing, and tail spars & ribs and 
fuselage bulkheads and walls 

S One variable 

• ^st_2nd-3rd-4th 

> Inner-wing spars & ribs 
•S One variable 

• 5 th 






% 


1 st : 0 ° 
2 nd : 45° 
3 rd : 90° 
4 th : -45° 


5 th : 0° 



6 th 

yth 


-45° 

90° 

45° 

0 ° 


■ 8 th 
gth 

Symmetric stacking of nine plies 
into a composite laminate 


Infeasible 



Structural components affected 
by thickness design variables 
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First Optimization Run 


Functions 

Performance indices 

Notes 

Objective 

FQQ = (PIw) 2 = W T 2 

DTOW 

Flutter 

constraint 

9j (X) = PI F = 1— x < 0- 

j = 1,2, ...,6 

15% margin 

Buckling 

constraint 

Qj (X) = PI B = (1/2) 2 — {positive min(BLF) — 1/2} 2 < 0. 

j = 7, 8,... ,11 

Safety factor =1.5 

Strength 

constraint 

g, (X) = PI S = — minp/IS) < 0. 
j = 12, 13,..., 16 

Safety factor =1.5 


□ Objective: total weight of gear up DTOW case 

□ Flutter constraints 

❖ Gear Up DTOW at M=0.66, 0.89, & 1.41 

❖ Gear up FFEP at M=0.66, 0.89, & 1.41 

□ Buckling & strength constraints 


❖ Minimum "buckling load factor" & minimum "margin of safety" from 
five analysis sets 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Load Cases 

1 

Up 

DTOW 

100, 200, 300, 600, 1100, 1200, 1300, 1400, & 1600 

2 

Up 

ZFW 

700, 800, 900, 1000, & 1700 

3 

Up 

M2W 

400 & 500 

4 

Down 

DTOW 

3001 -3017 + 3018 -3021 (emergency) + 1500 (for 
landing) 

5 

Down 

DLW 

4001 -4017 +4018 - 4021 (emergency) + 1800 (for 
landing) 
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O : Empty Fuel Empty Payload (EFEP) a : Full Fuel Empty Payload (FFEP) 
O: Zero Fuel Weight (ZLW) □: Design Take-Off Weight (DTOW) 


-10k S.L. 10k 20k 



First Optimization Run (continue) 



Performance Index 

Design Configuration 

Before Optimization 

Iteration 1 

Iteration 7 

Objective Function 

Total Weight 

DTOW; GU 

332738 

425764 

364105 

Constraint 

Functions 

gjOO 

Flutter 

Si 00 

DTOW; GU; M=0.66 

0.067(V) 

-0.362 

-0.342 

g 2 00 

DTOW; GU; M=0.89 

0.048fV) 

-0.543 

-0.096 

Ss 00 

DTOW; GU; M=1.41 

-0.079 

-1.34 

-0.297 

s + 00 

FFEP; GU; M=0.66 

0.066(V) 

-0.365 

-0.337 

Ss00 

FFEP; GU; M=0.89 

0.034(V) 

-0.586 

-0.094 

Se00 

FFEP; GU; M=1.41 

-0.095 

-1.32 

-0.255 

Buckling 

syOO 

DTOW; GU 

0.152fV) 

-1.05 

-1.29 

Sa00 

ZFW; GU 

0.186fV) 

-2.36 

-3.09 

S900 

M2W; GU 

0.151(V) 

-1.28 

-1.91 

Sio(X) 

DTOW; GD 

-0.960 

-3.27 

-3.88 

Sn00 

DLW; GD 

-0.561 

-0.308 

-1.07 

Strength 

S 12 O 0 

DTOW; GU 

0.999fV) 

-0.267 


s«00 

ZFW; GU 

0.998fV) 

-0.780 

f -0.061 \ 

SmOO 

M2W; GU 

0.997fV) 

-0.179 

l -0.537 J 

S 15 00 

DTOW; GD 

0.781fV) 

-0.751 

y5.63ey 

Si *00 

DLW; GD 

0.657fV) 

-0.210 


Weight penalty [%) 

0 

28.0 

9.43 1 


Manual 

update 


Design 

infeasible 


Design 

feasible 


Optimization based on DOT. 

Active and near active constraints 
from strength analyses. 
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First Optimization Run (continue) 


□ 



changes 


Mass balancing 
^ 2.0 


Upper skins 


Lower skins 


Spars, ribs, / 
bulkheads, & floors 


Fuselage skins 







|ll|f 

li 



Total thickness 
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Forward 


^ Color spectrums represent 
“thickness after”/“thickness before”. 
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First Optimization Run (continue) 



Active element 


□ Strain distribution of the active and near active constraints 




s 


Gear 


Forward 


down 


Main 


bay 


gear 


landing 



Main landing 
gear bay 



(b) Near active constraint (from strength 2; load case #1700) 


Main landing 
gear bay 


(a) Active constraint (from strength 4; load case #3013) (c) Near active constraint (from strength 1 ; load case #300) 
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Design Variables for the Second Optimization Run 




Forward 


Zone 2 Zone 3 

Zone 5 


Zone 6 


□ Composite ply angles of the six zones 

❖ Design variables are ply angles of the 
2 nd and 4 th layers. 

> Design variable linking 2 nd = -4 th 

❖ Use discrete design variables 


Bottom up 
view 


Zone 1 


Top down 

view Zone 1 


1 st : 0 ° 
2 nd : 45° 
3 rd : 90° 
4 th : -45° 


5 th : 0° 


6 th : -45° 
7 th : 90° 
8 th : 45° 
9 th : 0° 


Forward 


Value Ranges 

Discrete 

values 

0.0° < Ply angle < 2.5° 

0° 

2.5° < Ply angle < 7.5° 

5° 

7.5° < Ply angle < 12.5° 

10° 

12.5° < Ply angle < 17.5° 

15° 

17.5° < Ply angle < 22.5° 

20° 

22.5° < Ply angle < 27.5° 

25° 

27.5° < Ply angle < 32.5° 

30° 

32.5° < Ply angle < 37.5° 

35° 

37.5° < Ply angle < 42.5° 

40° 

42.5° < Ply angle < 47.5° 

45° 

47.5° < Ply angle < 52.5° 

50° 

52.5° < Ply angle < 57.5° 

55° 

57.5° < Ply angle < 62.5° 

o 

O 

VO 

62.5° < Ply angle < 67.5° 

65° 

67.5° < Ply angle < 72.5° 

70° 

72.5° < Ply angle < 77.5° 

75° 

77.5° < Ply angle < 82.5° 

o 

O 

oc 

82.5° < Ply angle < 87.5° 

85° 

87.5° < Ply angle < 90.0° 

90° 
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Second Optimization Run 


Functions 

Performance indices 

Notes 

Objective 

F(X) = — {0.5.<j, 2 (X) + 0.5 fll 3 (X) + 0 15 (X)} 

Safety factor =1.5 

Flutter 

constraint 

g, (X) = P, F = 1. - L1 g VL < 0. 
j = 1,2, ..., 6 

15% margin 

Buckling 

constraint 

gj (X) = PI B = (1/2) 2 — {positive min(BLF) — 1/2} 2 < 0. 

j = 7, 8,... ,11 

Safety factor =1.5 

Strength 

constraint 

gj (X) = PI S = -minpS) < 0. 
j = 12,13, 16 

Safety factor =1.5 


□ Purpose: Create more offset for active constraints 

□ Objective: performance index from the 2 nd and 4 th strength analysis sets 

□ Flutter constraints 


❖ Same as the first optimization run 
□ Buckling & strength constraints 


a DV, 


Active 
constraint 


Infeasible 

domain 



Tolerance 


Feasible 

domain 


Constraint 

boundary 


DVp Design Variable i 


♦♦♦ Same as the first optimization run except performance indices for the objective function 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Load Cases 

1 

Up 

DTOW 

100, 200, 300, 600, 1100, 1200, 1300, 1400, & 1600 

2 

Up 

ZFW 

700, 800, 900, 1000, & 1700 

3 

Up 

M2W 

400 & 500 

4 

Down 

DTOW 

3001 -3017 + 3018 - 3021 (emergency) + 1500 (for 
landing) 

5 

Down 

DLW 

4001 -4017 + 4018 - 4021 (emergency) + 1800 (for 
landing) 
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Second Optimization Run (continue) 



Performance Index 

Design Configuration 

Starting 

BBBC 1 

BBBC 2 

Objective Function 

-{ 

0.5 ai2 (X) + 0.5 5l3 (X) + ai5 (X)} 


0.111 

0.337 

0.348 



0i(X) 

DTOW; GU; M=0.66 

-0.342 

-0.341 

-0.342 


5_i 

<72 (X) 

DTOW; GU; M=0.89 

-0.096 

-0.096 

-0.096 


0) 
-M 
-1— > 

<73 (X) 

DTOW; GU; M=1.41 

-0.297 

-0.299 

-0.297 


Ft i 

<74® 

FFEP; GU; M=0.66 

-0.337 

-0.336 

-0.337 



05 (X) 

FFEP; GU; M=0.89 

-0.094 

-0.095 

-0.094 



06® 

FFEP; GU; M=1.41 

-0.255 

-0.257 

-0.255 

Constraint 

bJO 

<7? (X) 

DTOW; GU 

-1.29 

-1.50 

-1.38 

Functions 


<78 (X) 

ZFW; GU 

-3.09 

-3.63 

-3.09 

ofX) 

U 

<79 (X) 

M2W; GU 

-1.91 

-2.00 

-2.23 

y j vV 

3 

PQ 

<7io (X) 

DTOW; GD 

-3.88 

-4.22 

-4.19 



0n(X) 

DLW; GD 

-1.07 

-1.07 

-1.07 



012 (X) 

DTOW; GU 

-0.161 

-0.214 

-0.232 


-M 

bJO 

013 (X) 

ZFW; GU 

-0.061 

-0.141 

-0.145 


£ 

QJ 

014 (X) 

M2W; GU 

-0.537 

-0.204 

-0.542 


5-h 

-M 

GO 

015 (X) 

DTOW; GD 

-5.63e-6 

-0.159 

-0.159 



016 (X) 

DLW; GD 

-0.320 

-0.435 

-0.419 

Design Variables 

Starting 

BBBC 1 

BBBC 2 



1 (2 nd rib at inner-wing) 


45° 

15° 

15° 


2 (floor at main landing gear bay) 


45° 

65° 

65° 

3 (center wall at main landing gear bay) 

45° 

65° 

55° 

4 (aft bulkhead at main landing gear bay) 

45° 

65° 

55° 



5 (aft fuselage skin 1) 


45° 

30° 

40° 



6 (aft fuselage skin 2) 


45° 

50° 

55° 


Big-Bang Big-Crunch algorithm; number of population=60; number of Big-Bang Big-Crunch=2; discrete design variables 
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Third Optimization Run 



Third Optimization Run 


Functions 

Performance indices 

Notes 

Objective 

FQQ = (PIw) 2 = W T 2 

DTOW 

Flutter 

constraint 

9j (X) = PI F = 1— x < 0- 

j = 1,2, ...,6 

15% margin 

Buckling 

constraint 

Qj (X) = PI B = (1/2) 2 — {positive min(BLF) — 1/2} 2 < 0. 

j = 7, 8,... ,11 

Safety factor =1.5 

Strength 

constraint 

g, (X) = PI S = — minp/IS) < 0. 
j = 12, 13,..., 16 

Safety factor =1.5 


□ Objective: total weight of gear up DTOW case 


□ Flutter constraints 

❖ Gear Up DTOW at M=0.66, 0.89, & 1.41 

❖ Gear up FFEP at M=0.66, 0.89, & 1.41 

□ Buckling & strength constraints 


❖ Minimum "buckling load factor" & minimum "margin of safety" from 
five analysis sets 


Analysis 

Set 

Gear 

Configuration 

Weight 

Condition 

Load Cases 

1 

Up 

DTOW 

100, 200, 300, 600, 1100, 1200, 1300, 1400, & 1600 

2 

Up 

ZFW 

700, 800, 900, 1000, & 1700 

3 

Up 

M2W 

400 & 500 

4 

Down 

DTOW 

3001 -3017 + 3018 -3021 (emergency) + 1500 (for 
landing) 

5 

Down 

DLW 

4001 -4017 +4018 - 4021 (emergency) + 1800 (for 
landing) 
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Third Optimization Run 



Performance 

Index 

Design 

Configuration 

Starting 

Iteration 1 

Iteration 2 

Iteration 3 

Iteration 4 

Iteration 5 

Iteration 6 

Iteration 7 

Objective 

Function 

Total Weight 

DTOW; GU 

364105 

363810 







Constraint 

Functions 

gjOQ 

Flutter 

fli(X) 

DTOW; M=0.66 

-0.342 

-0.341 







flz(X) 

DTOW; M=0.89 

-0.096 

-0.096 







fl 3 (X) 

DTOW; M=1.41 

-0.297 

-0.297 







fl4(X) 

FFEP; M=0.66 

-0.337 

-0.336 







fls(X) 

FFEP; M=0.89 

-0.094 

-0.094 







.96 (X) 

FFEP; M=1.41 

-0.255 

-0.256 







Buckling 

.97 (X) 

DTOW; GU 

-1.38 

-1.38 







fl 8 (X) 

ZFW; GU 

-3.09 

-3.44 







.99 (X) 

M2W; GU 

-2.23 

-2.04 







.9io(X) 

DTOW; GD 

-4.19 

-3.95 







.9n(X) 

DLW; GD 

-1.07 

-1.07 







Strength 

.912 (X) 

DTOW; GU 

-0.232 

-0.227 







,9is(X) 

ZFW; GU 

-0.145 

-0.141 







,9i4(X) 

M2W; GU 

-0.542 

-0.261 







,9i 5 (X) 

DTOW; GD 

-0.159 

-6.16e-5 







■9l6(X) 

DLW; GD 

-0.419 

-0.359 







Weight penalty f %) 

9.43 

9.34 








► 


Optimization based on DOT. 
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Conclusions 


□ The Lockheed Martin's pre-matured N+2 LSCT aircraft is optimized in this study through the use of a multidisciplinary design 
optimization tool developed at the NASA AFRC. 

□ The baseline design of the pre-matured N+2 LSCT aircraft was infeasible when ZAERO based aeroelastic analyses were used. 

❖ This probably means that the aerodynamic loads distribution computed using ZAERO trim analysis are different than the MSC Nastran 
generated aerodynamic loads. 

□ The starting configuration of the optimization run should be an achievable design and weight penalty for this was 93,026 lb. 

❖ 28.0% increase from baseline 

□ During the first optimization run, the weight reduction was 61,659 lb, and therefore weight penalty at the end of the first 
optimization run is 31,367 lb. 

❖ Optimization was based on DOT optimizer 

❖ Active constraint: minimum margin of safety value is associated with the structural component located at the second rib of the inner wing near 
the main landing gear bay area. 

> Nose wheel yaw and steering case number 1 

❖ First near active constraint: minimum margin of safety value at the floor of main landing gear bay 

> 2.7g gust load case at Mach 0.89 and altitude of 20,000 ft 

❖ 9.4% increase from baseline 

❖ Second near active constraints: flutter speeds with DTOW and FFEP at Mach 0.89 

> Mass balancing effect to increase the flutter speeds 

□ The second optimization run was prepared to increase tolerance distance for the active and the first near active constraints. 

❖ Create more room for reducing total weight of the aircraft 

❖ Use six ply angles as design variables 

❖ Optimization was based on Big-Bang Big-Crunch algorithm with discrete design variables. 

❖ Can't change weight property, but can change strength property. Therefore, can create tolerance for future weight optimization run 
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Questions ? 




Future Studies 


□ Use lifting surface based aerodynamics 

❖ Use curvilinear sparibs to further reduce the weight of N+2 LSCT 

> Find % weight reduction through curvilinear sparibs technique 

❖ Add active control design variables 

> Use aeroelastic tailoring up to V L 

> Use active control between V L and 1.15V L 

> Find % weight reduction through game changing approach 

□ Use CFD based aerodynamics 

❖ Use more accurate air loads for optimizations 
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Backup: Object-Oriented MDO Tool 



Object-Oriented MDO tool 


❖ 


❖ 


❖ 


❖ 


❖ 


& 


Objective 
Function J &) 
Constraints 
G(x) 


□ Optimization is based on in-house Object-Oriented Optimization tool 

Equivalent to the following codes 

> Open MDAO, Model Center, Visual Doc, etc. 

Four optimizer codes are available. 

> Gradient based algorithms (Local optimizers) 

S DOT 
y ADS 

> Global optimizers (Gradient free algorithms) 

y Genetic Algorithm 
y Big-Bang Big-Crunch Algorithm 

□ Update design pre-processor module 

❖ Update MSC/N ASTRAN input file 

□ Modal analysis module 
Perform modal analysis using MSC/NASTRAN sol. 103 
Save following data 

> Total weight, CG location, mass moment of inertia 

> Frequencies & mode shapes and global mass matrix 

□ Weight post-processor module 

Use MSC/NASTRAN sol 103 results for small weight. 

> MSC/NASTRAN results has number of digit issue. 

Use in-house weight computation code for large weight. 


Optimizer 


Object 
Oriented 
Optimization 
Tool 


Si 


Design 

Variables 






proce 

#i 

1 

icipl 

#i 



Need to develop; 
Depends on 
Analysis code 


processor p ost _ 
"J nronessi 


processor 


Discipline 
#j 



[Performance] 

/ Script ] / 

A Indices / 

l Commands// | 


mb 


Pre- 
processor 
#k Post- 
Discipline processor 
#k 


Analysis code; 
Commercial and/or 
in-house codes 



❖ 
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Object-Oriented MDO tool (continue) 


□ Flutter analysis and flutter post-processor modules 

❖ Use ZAERO code for flutter analyses 

❖ Use an in-house flutter speed tracking program 



V F > 1.15V l 




□ Update ZAERO pre-processor, trim analysis, trim loads pre-processor modules 

❖ Update ZAERO input data 


> Based on total weight, CG locations, moment of inertias, and global mass matrix 

❖ Use ZAERO code for trim analysis 

> Create design loads for various design configurations 

❖ Post-process the splined loads 

> Create symmetric and anti-symmetric loads 


□ Landing and ground control loads pre-processor module 

❖ Compute corresponding design loads using in-house code 

> Landing loads 

> Ground control loads 

> Emergency landing loads 
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Object-Oriented MDO tool (continue) 

□ Buckling and strength analyses and strength post-processor modules 

❖ Based on MSC/NASTRAN sol. 105 

❖ Use in-house strength post-processor code 

❖ Safety factor of 1.5 is used for all metal and composite materials in this study. 


Design Load x Safety Factor < Failure Load 


Failure Load 

Design Load x Safety Factor 


< 0 


MS = 


Failure Load 

Design Load x Safety Factor 


□ Buckling post-processor module 

❖ Use in-house code 

❖ Buckling Load Factor (BLF) 0<BLF<1 : Buckling predicted BLF<0 or BLF>1: Buckling not predicted 

> Buckling predicted: 0 < BLF < 1 - 1/2 < BLF - 1/2 < 1/2 (BLF - 1/2) 2 < (1/2) 2 


> Buckling not predicted: (BLF — 1/2) 2 > (1/2) 2 


(1/2) 2 — ( BLF — - 


< 0 . 


PI B = (1/2) 2 — (positive min(BLF) — 1/2} 2 


PI S = -min(MS) 
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